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. ABSTRACT 

o 

The stellar upper-mass limit is highly uncertain. Some studies have claimed there is a universal upper limit of ~150Mg. A factor 
^..^1 that is often overlooked is that there might be a significant difference between the present-day and the initial masses of the most 

. massive stars - as a result of mass loss. The upper-mass limit may easily supersede ~200Mo. For these reasons, we present new 

mass-loss predictions from Monte Carlo radiative transfer models for veiy massive stars (VMS) in the mass range 40-300 Mq, and 
with very high luminosities 6.0 < log(L*/LQ) < 7.03, corresponding to large Eddington factors F. Using our new dynamical approach, 
we find an upturn or "kink" in the mass-loss versus F dependence, at the point where the model winds become optically thick. This 
coincides with the location where our wind efficiency numbers surpass the single-scattering limit of = 1, reaching values up to 
^ 2.5. In all, our modelling suggests a transition from common 0-type winds to Wolf-Rayet characteristics at the point where the 
winds become optically thick. This transitional behaviour is also revealed with respect to the wind acceleration parameter, yS, which 
I starts at values below 1 for the optically thin O-stars, and naturally reaches values as high as 1.5-2 for the optically thick Wolf-Rayet 

models. An additional finding concerns the transition in spectral morphology of the Of and WN characteristic He n line at 4686A. 
When we express our mass-loss predictions as a function of the electron scattering Eddington factor Fc ~ L^/Mt, alone, we obtain 
' an M vs. Fe dependence that is consistent with a previously reported power law M oc F^ (Vink 2006) that was based on our previous 

CIh, semi-empirical modelling approach. When we express M in terms of both Fe and stellar mass, we find optically thin winds and M 

oc M^^'^^rl'- for the Fe range 0.4 :S Fe S 0.7, and mass-loss rates that agree with the standard Vink et al. recipe for normal O stars. 
For higher Fe values, the winds are optically thick and, as pointed out, the dependence is much steeper, M oc Mt'^''^rf''^ . Finally, 
we confirm that the effect of F on the predicted mass-loss rates is much stronger than for the increased helium abundance (cf. Vink 
^ . & de Koter 2002 for Luminous Blue Variables), calling for a fundamental revision in the way stellar mass loss is incorporated in 

evolutionary models for the most massive stars. 
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> 

' 1. Introduction rates were obtained that were a factor ~3 higher than for cases 

in which single scattering was strictly enforced. 

O The prime aim of this paper is to investigate the mass-loss be- Historically, the situation for the WR stars was even more ex- 

• . haviour of very massive stars (VMS) with masses up to 300 treme. Here, t, values of -10 had beenreported (e.g. Barlow et al. 

^ ; Mo that approach the Eddington limit. Mass loss from hot mas- ^g^^y ^^jj^ identification of major wind-clumping eff^ects on 

sive stars is driven by radiative forces on spectral lines (Lucy & ^j^^ ei^pirical mass-loss rates (Hillier 1991; Moffat et al. 1994; 

Solomon 1970; Castor, Abbott & Klein 1975; CAK). CAK de- Hamann & Koesterke 1998), these numbers should probably be 

veloped the so-called force multiplier formalism in order to treat down-revised to values of 77 -3. Although it has been argued 

all relevant ionic transitions. This enabled them to simultane- ^j^^^ ^-^^^ ^j^^ ^^^^^ radiation pressure (Lucy & 

■ ously predict the wind mass-loss rate, M, and termmal velocity, ^^bott 1993, Springmann 1994, Gayley & Owocki 1995, Nugis 
'^C ; Of O-type stars. Although these predictions provided reason- ^ tamers 2002, Grafener & Hamann 2005), the prevaiUng no- 
_ able agreement with observations, they could account for neither ^^^^ ^^^jj ^^^^ ^^^^^ ^pjj^^jjy ^^^^^ ^^^^^^^ ^^^^^ ^^^^^ 

the high wmd efficiencies r] = Mv^/L/c of the denser Of stars ^^^-^ ^^-^^ accelerating flow lies within the pseudo or 

with their strong He 11 4686A lines, nor that of the even more ^ise photosphere, are fundamentally different from the transpar- 

extreme Wolf-Rayet (WR) stai's. ent line-driven O-star winds (e.g. Grafener & Hamann 2008) . 

This discrepancy had been proposed as due to the neglect of por O-type stars Miiller & Vink (2008) recendy suggested 
multi-line scattering (Lamers & Leitherer 1993, Puis et al. 1996). a new parametrization of the line acceleration, expressing it as a 
Using a global energy Monte Cai'lo approach (Abbott & Lucy function of radius rather than of the velocity gradient (as in CAK 
1985, de Koter et al. 1997) in which the velocity law was adopted theory). The implementation of this new formalism improves the 
- aided by empirical constraints - Abbott & Lucy (1985) and local dynamical consistency of Monte Carlo models that origi- 
Vink et al. (2000) provided mass-loss predictions for galactic O nally imposed a velocity law. Not only do we find fairly good 
stars including multi-Une scattering. This appeared to solve the agreement with observed terminal velocities (see also Muijres et 
windmomentumproblemforthedenserO-starwinds. Mass-loss al. 2011b), but as our method naturally accounts for multi-line 
scattering, it is also applicable to denser winds, such as those of 
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Still adopting a velocity stratification, Vink & de Koter 
(2002) and Smith et al. (2004) predicted mass-loss rates for 
Luminous Blue Variables (LBVs), and showed that M is a strong 
function of the Eddington factor for these objects. They also 
found that, despite their extremely large radii, even LBV winds 
may develop pseudo-photospheres under special circumstances: 
when they find themselves in close proximity to the bi-stability 
and Eddington limits. 

In this paper, our aim is to study the mass-loss behaviour 
of stars as they approach the Eddington limit (see also Grafener 
& Hamann 2008). We do this in a systematic way by targeting 
VMS in the range 40-300 Mq. A pilot study was performed by 
Vink (2006) who found a steep dependence of M on Fg, finding 
M oc Eg, but this was obtained using the earlier global energy 
approach, in which the velocity stratification was adopted, rather 
than our new dynamically-consistent approach that we explore 
in the following. 

Both approaches have their pros and cons. In the semi- 
empirical approach, the terminal velocity constraint (where 
yoo/uesc is constant) might aid the modelling at imposing the cor- 
rect wind structure, and as long as the adopted velocity law is 
close to the correct one, it provides the most accurate mass-loss 
rates. The prime advantage of this approach is that any miss- 
ing physics that could affect the mass-loss rate might be bal- 
anced by employing the empirical (and probably close to real- 
istic) terminal-wind velocity. The second approach, however, is 
theoretically more appealing, as it enforces local dynamical con- 
sistency, and one thus no longer needs to rely on free parameters. 
Ultimately, one would aspire to predict accurate mass-loss rates 
and terminal velocities simultaneously from first principles via 
the second approach explored here. 

The stellar upper-mass limit is highly controversial. On 
purely statistical grounds, some investigators have claimed there 
is a universal upper limit of ~15OM0 (e.g. Weidner & Kroupa 
2004, Oey & Clarke 2005, Eiger 2005). However, a physical fac- 
tor that is often overlooked concerns the possibility of a signifi- 
cant difference between the present-day and the initial masses of 
the most massive stars, as a result of strong mass loss. In other 
words, the initial masses of the most massive stars may be signif- 
icantly above I5OM0, possibly superseding 200 Mq (e.g. Eiger 
et al. 1998, Crowther et al. 2010, Bestenlehner et al. 201 1). The 
issue of the upper mass-limit will remain uncertain as long as 
there is only limited quantitative knowledge of mass loss in close 
proximity to the Eddington limit. 

Our aim is thus to explore wind models of stars with masses 
up to 3OOM0, using a well-established methodology that has 
been extensively tested against observations for lower mass 
common O-type stars. VMS have been proposed as leading to 
the production of intermediate mass (of the order of IOOMq) 
black holes that have been suggested to be at the heart of ultra- 
luminous X-ray sources (Belkus et al. 2007 and Yungelson et al. 
2008). Clearly, the success of such theories depends critically on 
the applied mass-loss rates. The present study may help advance 
these theories. 

Our paper is organized as follows. In Sect. |2]we briefly de- 
scribe the Monte Carlo mass-loss models, before presenting the 
parameter space considered in this study (Sect.|3). The mass-loss 
predictions (Sect. |4]i are followed by a description of the spec- 
tral morphology of the Of-WN transition in terms of the char- 
acteristic He II 4686A line in Sect. |5] Subsequently, we compare 
our wind model parameters against empirical values for the most 
massive stars in the Arches cluster (Sect. lSTt . as well as theoret- 
ical models (Sect.jOJ of CAK and Grafener & Hamann (2008), 
before ending with a summary in Sect. |7] 



2. Monte Carlo models 

Mass-loss rates are calculated with a Monte Carlo method that 
follows the fate of a large number of photon packets from be- 
low the stellar photosphere throughout the wind. The core of 
our approach is related to the total loss of radiative energy that 
is coupled to the momentum gain of the outflowing material. 
Since the absorptions and scatterings of photons in the wind de- 
pend on the density in the wind, hence on the mass-loss rate, 
one can obtain a consistent model where the momentum of the 
wind material equals the transferred radiative momentum. We 
have recently improved our dynamical approach (Miiller & Vink 
2008, Muijres et al. 2011b) and are now able to predict M si- 
multaneously with Vca and the wind structure parameter /?. The 
essential ingredients and assumptions of our approach have been 
discussed more extensively in Abbott & Lucy (1985), de Koter 
et al. (1997), and Vink et al. (1999). Here we provide a brief 
summary. 

The Monte Carlo code mc-wind uses the density and tem- 
perature stratification from a prior model atmosphere calcula- 
tion performed with isa-wind (de Koter et al. 1993, 1997). These 
model atmospheres account for a continuity between the pho- 
tosphere and the stellar wind, and describe the radiative trans- 
fer in spectral lines by adopting an improved Sobolev treatment. 
The chemical species that are explicitly calculated (in non-LTE) 
are H, He, C, N, O, S, and Si. The iron-group elements, which 
are crucial for the radiative driving and the M calculations, are 
treated in a generalized version of the "modified nebular approx- 
imation" (e.g. Schmutz 1991). However, we performed a number 
of test calculations in which we treated Ee explicitly in non-LTE. 
These tests showed that differences with respect to the assump- 
tion of the modified nebular approximation for Ee were small. 
Therefore, we decided to treat Ee in the approximate way, as 
was done in our previous studies. 

The line list used for the MC calculations consists of over 10^ 
of the strongest transitions of the elements H - Zn extracted from 
the line list constructed by Kurucz & Bell (1995). The wind was 
divided into 90 concentric shells, with many narrow shells in the 
subsonic region, and wider shells in supersonic layers. Eor each 
set of model parameters, a certain number of photon packets was 
followed, typically 2 lO*". 

Other assumptions in our modelling involve wind stationar- 
ity and spherical geometry. The latter seems to be a good approx- 
imation, as the vast majority of O-type stars show little evidence 
of significant amounts of linear polarization (Harries et al. 2002, 
Vink et al. 2009). Nevertheless, asphericity has been found in 
roughly half the population of LBVs (Davies et al. 2005, 2007), 
although those polarimetry results have been interpreted as the 
result of small-scale structure or "clumping" of the wind, rather 
than of significant wind asymmetry. 

With respect to wind clumping, it has been well-established 
that small-scale clumping of the outflowing gas can have a pro- 
nounced effect on the ionization structure of both O-star and 
Wolf-Rayet atmospheres (e.g. Hillier 1991). This has lead to a 
downward adjustment of empirical mass-loss rates, by factors of 
up to three (e.g. Moffat et al. 1994, Hamann & Koesterke 1998, 
Mokiem et al. 2007, Puis et al. 2008), and possibly even more 
(Bouret et al. 2003; Eullerton et al. 2006). In addition, clump- 
ing may have a direct effect on the radiative driving, therefore 
on predicted mass-loss rates (e.g. Grafener & Hamann 2008). 
The subtle issues of both clumping and porosity on the predicted 
mass-loss rates have recently been investigated by Muijres et al. 
(201 la). Whilst it was found that the impact on M can be high for 
certain clumping prescriptions, the overall conclusion was that 
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clumping does not affect the wind properties of O-type dwarfs 
and supergiants in a dramatic way for moderate clumping fac- 
tors and porosity. Stars close to the Eddington limit, however, 
may be much more susceptible to - even modest - clumping 
(Shaviv, 1998, 2000; van Marie et al. 2008). Another relevant 
factor that might affect the results concerns the interplay of ra- 
diation with rotation, as the Eddington factor for a rotating star 
depends explicitly on the rotation rate (Langer 1997, Maeder & 
Meynet 2000). As a consequence, the effective Eddington limit 
could be reduced by rotation, and might even become the domi- 
nant factor 

In the present set of computations, we do not account for 
the effects of wind clumping or rotation, but it should be kept in 
mind that these effects might play a quantitative role. 

3. Parameter space and model applicability 

Stars approach the Eddington limit when gravity is counterbal- 
anced by the radiative forces; i.e. F = grad/^Newton = 1- Photons 
can exert radiative pressure through bound-free, free-free, elec- 
tron scattering, and bound-bound interactions. In early-type stars 
hydrogen, the dominant supplier of free electrons, is fully ion- 
ized. Therefore Fe - /^Newton is essentially independent of 
distance and constitutes a fixed number for each model. Because 
of this useful property, which provides a well-defined and simple 
quantitative handle, we opt to discuss our results in terms of Fe. 
We discuss this choice in more detail in Sect. 13.11 

The dependence of the mass-loss rate M on Fe represents a 
non-trivial matter because M depends on both the mass M and 
the stellar luminosity L. To properly investigate the effect of high 
Fe on mass-loss predictions, we first need to establish the rele- 
vant part of parameter space in terms of M, L, and Fg. For a fully 
ionizecQ plasma, Fg equals 

re = ^^ = -^ = 10-^-«'^(l+X)(^(^r , (1) 

The luminosities are chosen in such a way that in combination 
with the stellar mass M the desired Fe value is obtained. It is 
principally the effective temperature that sets the ionization strat- 
ification in the atmosphere, and thus determines which lines are 
most active in driving the wind. As a result, T^g affects the pre- 
dicted mass-loss rate. For most parts of this paper, we investigate 
the influence of Fe for a fixed stellar temperature of 50 000 K. 
The Teff dependence is studied separately in Sect. 14.41 All mod- 
els are for the solar metallicity from Anders & Grevesse (1989) 
and with the element-to-element distribution from Allen (1973). 
The prime reason we use these older abundances rather than the 
newer (and lower) solar or B-type star abundances is to be able 
to directly compare the new results to the older Vink et al. (2000) 
rates. We note however that it is the element of Fe that sets 
the mass-loss rate, and this element has not changed. The low 
(CNO) and intermediate-mass elements however dominate the 
outer wind, where the terminal wind velocity is set. Nonetheless, 
even a substantial decrease in these abundances (by several tens 
of percent) is not expected to lower the terminal wind velocities 
significantly, since the terminal wind velocity has been found to 
depend only weakly on metallicity (Leitherer et al. 1992). 

We divided our model stars into three different groups ac- 
cording to their characteristics. The first group comprises ob- 
jects that have relatively common O-star masses in the range 

' In reality, Fe changes once hydrogen recombines (which starts be- 
low 30000 K), or when the hydrogen-to-helium surface abundance 
changes, relevant for classical WR stars. 



40-70 Mq. The second group of objects are rather high-mass 
stars within the "observable range" of 70-120 Mq. They might 
be close to the Eddington limit already early-on on the main se- 
quence because of their intrinsically high luminosity. The third 
group involves very massive stars in the mass range 120-300 
Mq. They are near the Eddington limit for the same reason as 
the second group. So far, there is a lack of compelling obser- 
vational evidence of any such stars in the present-day universe; 
however, we note that Crowther et al. (2010) have suggested a 
revision of the upper-mass limit to ~300 Mq. 

The bulk of the models in our grid have been chosen such 
that the behaviour of mass loss as a function of M and L can 
be studied separately. The grid is presented in Table [1] We note 
that the (M, L) combinations are intentionally rather extreme to 
assure high Fe values. The reason is to specifically map that part 
of parameter space where physically the most extreme winds are 
expected to appear. 

3. 1 . Model applicability regime 

With respect to the potential limitations of our modelling ap- 
proach, we make one rather stringent assumption in the manner 
the (sub-) photospheric density structure is set-up. In the deepest 
layers of the model atmosphere (with v << Ikm/s), we assume 
that the run of density is provided by the equation of motion us- 
ing grad - ge, SO wc apply F = Fe. In reality F > Fe, as well as 
being depth-dependent as a result of bound-bound, bound-free, 
and free-free processes. Notably, the opacities from millions of 
weak iron lines may contribute significantly, but they are largely 
neglected in the deep layers of our models. 

Nugis & Lamers (2002) highlight the importance of the iron 
peak opacities in deep photospheric layers for the initiation of 
Wolf-Rayet winds (see also Heger & Langer 1996). This ap- 
proach was subsequently included in models by Grafener & 
Hamann (2005, 2008) for WC and WNL stars. They find that 
the presence of these opacity bumps may locally cause F to ap- 
proach unity, leading to the formation of optically thick winds. 
In our Monte Carlo approach, we trace the radiative driving of 
the entire wind, and as most of the energy is transferred in the 
supersonic part of the outflow, we are less susceptible to the de- 
tails of the (sub)photospheric region. However, this also means 
that we do not treat these deep regions self-consistently. This 
implies that we can (and we will) compute model atmospheres 
with values of Fe very close to one. This strategy has the advan- 
tage of allowing us to explore the transition from transparent to 
dense stellar winds. As our models do capture the full physics in 
the layers around and above the sonic point, we argue that they 
correctly predict the qualitative behaviour of dense winds, but 
that Fe for one of our optically thick wind models would corre- 
spond to a model with smaller Fe if the ionic contributions were 
included in the deepest parts of the atmosphere. This "shift" in 
Fe is not fixed but would depend on the sonic point temperature 
and density. From the behaviour of the Rosseland mean opacity, 
we would expect the size of the shift to increase at higher Fe and 
higher temperatures. 

If F exceeds unity at some depth in the subphotospheric part 
of the atmosphere, a density inversion is expected to occur for 
the static case, i.e. for increasing radial distance from the centre, 
the density very near the domain where F > 1 is anticipated to in- 
crease. This is encountered in studies of stellar structure and evo- 
lution, but it is unclear what really happens in nature. The poten- 
tial effects may involve strange-mode pulsations (e.g. Glatzel & 
Kiriakidis 1993), subsurface convection (Cantiello et al. 2009), 
or an inflation of the outer stellar envelope (e.g. Ishii et al. 1999). 
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Fig. 1. The predicted mass-loss rates versus Tg for models ap- 
proaching the Eddington limit. Asterisks, diamonds, and trian- 
gles correspond to models of the respective mass ranges I, II, 
and III. Our model assumptions likely break down to the right of 
the vertical dashed Une. 



Fig. 2. The logarithmic difference between the new Fg mass-loss 
predictions and the standard Vink et al. (2000) recipe for models 
approaching the Eddington limit. Symbols are the same as in 
FigH 



These processes tend to occur only when F is very close to unity 
or above it (see e.g. Petrovic et al. 2006). In assessing the out- 
come of our computations, we find that at Fg > 0.95 the results 
behave rather oddly. Though we show the wind results over the 
entire Fg range, we only quantify the mass-loss rates up to this 
value of Fg. This boundary is indicated by a vertical dashed line 
in all relevant figures. 



4. Results 

4.1. Mass-loss predictions at higli Fg 

Table \T\ lists our mass-loss predictions for all three considered 
mass ranges. Most columns are self-explanatory, but we note 
that the effective escape velocity fgsc'^'* (7th column) is defined as 
V2GM(1 - Fe)/7?. The predicted wind terminal velocities, mass- 
loss rates, wind efficiency numbers, and wind acceleration pa- 
rameter jS are given in columns (8), (9), (10), and (11), respec- 
tively. For comparison column (12) lists the mass-loss values 
from the standard mass-loss recipe of Vink et al. (2000) where /3 
was held fixed at unity. 

The predicted mass-loss rates (column 9) are shown in Fig.[T] 
Dififerent symbols are used to identify the different mass ranges. 
The figure shows that M increases with Fg. This is in qualitative 
agreement with the luminosity dependence of the standard mass- 
loss recipe of Vink et al. (2000), derived from a set of models 
with Fg ^ 0.4. Analogous to the results from the standard Vink et 
al. (2000) recipe, Fig.[T]suggests that there is an additional mass- 
loss dependence on mass, as for fixed Fg the higher mass stars 
have higher mass-loss rates. This finding confirms that mass-loss 
rates cannot solely be described by a dependence on luminosity 
or Eddington factor. This will be discussed further in Sect. 14.21 

When comparing columns (9) and (12) from Table [1] it can 
be noted that our new high Fg mass-loss predictions tend to 
be larger than those determined using the standard Vink et al. 
(2000) recipe. In order to quantify these differences, we divide 
the new mass-loss rates over those determined using the Vink et 
al. (2000) recipe (using the derived terminal wind velocities as 
input), and show the results in Fig. |2] For the range Fg $0.7, the 
differences are small. However, for values of Fg exceeding ~0.7, 
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Fig. 3. The predicted terminal wind velocities versus Fg for mod- 
els approaching the Eddington limit. Symbols are the same as in 
FigH] 



the new and the old results diverge sharply. The maximum dif- 
ference reaches a factor of five, which is similar in magnitude to 
what was reported previously for LBVs (Vink & de Koter 2002) 
and WR stars (Vink & de Koter 2005). We note that, although 
these prior results were based on global energy consistency with 
fixed Woo/ygsc, where the velocity stratification was adopted, the 
reason for the differences revealed in Fig. |2]is that we probe a 
different part of parameter space. 

We now turn our attention to the wind velocity structure. We 
first inspect the associated terminal wind velocity predictions. 
Figure [3] shows the behaviour of terminal wind velocity versus 
Fg. The highest values are reached for the highest mass stars 
and exceed 5000 km/s. As expected drops with Fg. In the 
Fg range 0.4-0.95, the terminal wind velocity divided over the 
escape velocity is of the order 3-4, which is similar to the values 
for common O-type stars (Muijres et al. 201 lb), where T^s is in 
the range 30-40 kK, and the wind velocities are closer to 3000 
km/s (see Sect. [441. 
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1240 


351 


3595 


-4.33 


2.17 


1.60 


-4.74 


20 


90 


6.50 


0.93 


23.7 


1203 


318 


3613 


-4.40 


2.29 


1.60 


-4.81 


21 


100 


6.57 


0.982 


25.7 


1218 


163 


(1752) 


(-3.81) 


(2.78) 


(2.10) 


-4.46 


22 


85 


6.50 


0.987 


23.7 


1169 


133 


(1808) 


(-3.89) 


(3.11) 


(2.09) 


-4.57 


Mass Range III: 
























23 


120 


6.42 


0.58 


21.6 


1455 


945 


5744 


-5.00 


1.03 


1.17 


-4.96 


24 


120 


6.50 


0.70 


23.7 


1389 


761 


5122 


-4.78 


1.29 


1.22 


-4.82 


25 


300 


6.97 


0.83 


40.7 


1676 


700 


5527 


-4.20 


1.74 


1.34 


-4.36 


26 


180 


6.76 


0.85 


31.95 


1465 


568 


4642 


-4.31 


1.76 


1.37 


-4.53 


27 


250 


6.91 


0.86 


38.0 


1583 


593 


4885 


-4.16 


1.87 


1.37 


-4.40 


28 


225 


6.87 


0.87 


36.3 


1537 


554 


4657 


-4.17 


2.03 


1.43 


-4.43 
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1 .HO 
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30 


300 


im 


0.95 


43.6 


1619 


362 


(3728) 


(-3.82) 


(2.34) 


(1.62) 


-4.23 


31 


200 


6.87 


0.977 


36.3 


1449 


219 


(2500) 


(-3.63) 


(2.51) 


(2.17) 


-4.30 


32 


120 


6.65 


0.986 


28.1 


1276 


151 


(3136) 


(-4.10) 


(2.87) 


(1.68) 


-4.69 


33 


250 


6.96 


0.97 


40.2 


1539 


267 


(4085) 


(-3.96) 


(2.82) 


(1.93) 


-4.40 


34 


153 


6.75 


0.975 


31.6 


1359 


215 


(3409) 


(-4.02) 


(3.04) 


(1.81) 


-4.58 



Table 1. High Fe mass-loss predictions for all three mass-range grids, with T^s kept constant at 50,000 K. The stellar parameters 
are given in Columns (2-7), providing the stellar mass, luminosity, Eddington factor, radius, escape velocity, and the effective 
escape velocity. Columns (8-11) give the wind properties: the terminal velocity, the mass-loss rate, the wind-efficiency number 
J] — MvcolLi,c, and the wind acceleration parameter The last column provides the mass-loss rates as computed using the formula 
by Vink et al. (2000) with the computed terminal wind velocity as input. 



We next turn our attention to the other wind velocity struc- 
ture parameter, yS, which describes how rapidly the wind accel- 
erates. The predicted values of fi are depicted in Fig. |4] but fi 
does not show a significant dependence on stellar mass. For Fg 
up to 0.7, fi values are near unity, in accordance with the dy- 
namical consistent models of Pauldrach et al. (1986), Miiller & 
Vink (2008), and Muijres et al. (2011b). However, when Fg ex- 
ceeds 0.7 and approaches unity,/? steadily rises to values of about 
1.7. These higher p values are supposedly more commensurate 
in Wolf-Rayet stars (see e.g. Ignace et al. 2003), and it is reas- 
suring to find that our models naturally predict this transition, 
without the use of any free parameter. 

In all, our results suggest a natural extension from 0-type 
mass loss to more extreme WR behaviour for increasing Fg. An 
upturn in the M behaviour is found at Fg ~0.7. Inspection of our 
models reveals a change from optically thin to optically thick 
wind models at the position where we obtained the kink in the 
mass-loss versus F relationship. Closer scrutiny of our model 



output also revealed that here the character of the Fe line driving 
changes. Whilst various ionization stages of Fe contribute at all 
Fg models, we find that for the optically thin models at the low 
Fg end, just one ionization state of Fe dominates the relevant part 
of the wind driving domain (from just below the sonic point to 
about half the terminal velocity). By contrast, for the optically 
thick models at the high Fg end, two or more ionization stages of 
Fe contribute to the primary driving regime. 

4.2. Fg dependence of mass loss 

In order to determine the dependence of the mass-loss rate on 
Fg, we could simply fit the datapoints to a power law: 

M oc Fg'' (2) 

Using the semi-empirical approach Vink (2006) found p to be 
equal to ~5, and our dynamically consistent results provide the 
same slope here. However, in order to also take the mass de- 
pendence into account we divide the mass-loss rates by M^', and 
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Fig. 4. The predicted wind velocity structure parameter /? versus 
Te for models approaching the Eddington limit. Symbols are the 
same as in FigH] 
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Fig. 5. The predicted mass-loss rates divided by versus Tg 
for models approaching the Eddington limit. The dashed-dotted 
line represents the best linear fit for the range 0.4 < 0.7. The 
dashed line represents the higher 0.7 < 0.95 range. Symbols 
are the same as in FiglT] 



M oc M^-^'^Y^-'^ 



For 0.7 < rp< 0.95 we determine that 



0.78t- 4.77 



MocM' 

We first note that we intentionally do not provide equations 
here, as we expect the absolute values of these mass-loss predic- 
tions for the high T models to be underpredicted. In Sect. 16.11 
we see that our predicted wind terminal velocities are a factor 
2-4 higher than found empirically. If we had employed our older 
semi-empirical approach, and assumed empirical - i.e. a factor 
2-4 lower - terminal wind velocities, we would have predicted 
higher mass-loss rates. Indeed, the pilot study of Vink (2006) 
that was based on the semi-empirical approach provided such 
higher mass-loss rates. 

Secondly, we note that the exact value of the Fg transition 
value is model-dependent. We have already discussed our po- 
tential model deficiencies in the deepest layers in Sect. 13.11 but 
we should note that several other factors might also play a role; 
in particular, higher mass-loss rates - as would be obtained from 
our semi-empirical approach - would shift the kink to lower Fg 
values. We note that the latter could occur in the case the wind is 
clumped. As long as porosity effects are small, wind clumping 
is expected to increase the momentum transfer in our MC mod- 
els (Muijres et al. 2011a). Although future dynamical-consistent 
modelling of clumped winds is required to test this, wind clump- 
ing could potentially increase our predicted mass-loss rates, and 
subsequently decrease the Fg value of the kink. 

The above mass-loss relationships can easily be transformed 
using Eq. ([T]i, leading to the entirely analogous mass-loss re- 
lationships M cc LO-6'^Fe'-^^ and Eq. © to M cc L^-^'^Fe^ ''^ 
Interestingly, if one subsequently applies a mass-luminosity re- 
lationship for classical (He-rich) WR stars of Maeder & Meynet 
(1987) or for very massive H-rich stars such as that of Yungelson 
et al. (2008), with L oc M'^"* for both cases, it follows that 
M oc M^ "*. This appears to be in good accord with the radio 
mass-loss rate relation M oc M^^ for classical WR stars with 
measured masses from binaries by Abbott et al. (1986). It also 
agrees with the M versus stellar mass relationship of M oc M"^ 
that has been applied in WR evolution models by Langer (1989). 



show the results in Fig. |5] We fit the data with the following 
power-law 

M oc M^T/ (3) 

Below Fe ^ 0.95, Fig. |5] shows two mass-loss regimes, divided 
by a boundary at Fe ~ 0.7. This is not only the point where the 
slope of the mass-loss versus F relation changes, but also where 
the wind efficiency parameter 77 surpasses the single scattering 
limit (see below). Upon further inspection of our models, we find 
that as long as Fe ^ 0.7 the winds are optically thin, implying 
that the sonic point of the outflowing material lies outside the 
photosphere, whilst the winds become optically thick - with the 
photosphere moving outside of the sonic point - for Fg values 
above ^0.7. 

We derive two independent mass-loss relationships for the 
two separate Fe regimes. 

For0.4<Fe<0.7 we find 



4.3. Increased wind efficiency close to the Eddington limit? 

In order to learn whether radiation-driven mass-loss rates con- 
tinue to increase with increasing Fe or reach a maximum in M 
instead, it is insightful to consider the wind efficiency parameter 
77 - MvooliLIc). We show the predicted values of 77 in Fig.|6] As 
the symbols denote different mass ranges, the small scatter on 
the datapoints shows that i] is not very sensitive to stellar mass. 
At values of Fe ~ 0.5 we find wind efficiency numbers 77 of or- 
der 1, in accordance with standard Vink et al. (2000) models. 
However, when Fe approaches unity, 77 rises in a curved man- 
ner to values as high as 77 ^ 2.5. Such high 77 values are more 
commensurate with Wolf-Rayet winds than with common O star 
winds, and these results thus confirm a natural extension from 
common O-type mass loss to more extreme WR behaviour. In 
Sect. 16. II we find that our predicted wind-terminal velocities are 
higher than the empirical values. As we note that an overpre- 
diction of the wind velocity is likely offset by a mass-loss rate 
underprediction by a similar amount, we argue that the combi- 
nation of these two quantities, i.e. their product constituting the 
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Fig. 6. The predicted wind efficiency number rj versus Te for 
models approaching the Eddington limit. Symbols are the same 
as in FiglT] 



T] parameter, might be less affected by model deficiencies than 
either of these quantities would be individually. 

The maximum mass loss in our models up to Fg - 0.95 is 
log Mm-dx - -3.8. This is the mass-loss rate that is retrieved for 
the most extreme models in our grid. Owocki et al. (2004) have 
investigated the mass loss of stars that formally exceed their 
Eddington limit and show that the expected mass loss falls well 
below the values required to account for the mass that is lost 
during LBV giant eruptions, such as that of t] Carinae in the 
1840s. Interestingly, they introduce a porosity-moderated con- 
tinuum driven mass loss that might account for the huge mass- 
loss rates associated with LBV eruptions (which may be of the 
order of IMo/yr). 



Fig. 7. The predicted mass-loss rates versus effective tempera- 
tures for several values of Tg, with from top to bottom Fg equal 
to 0.90 (model 29; open square), 0.87 (model 28;open triangle), 
0.83 (model 25;open diamond), and 0.58 (model 23; asterisk), 
respectively. 



4.4. Effect of Teff on high Fg models 

To establish whether there is an additional temperature depen- 
dence on M, we varied Tgff over the range 50-30 kK for selected 
Fg models, with mass-loss predictions presented in Fig. |7] and 
terminal wind velocities shown in Fig. |8] As we wish to stay 
above the temperature of the bi-stability jump (which starts at 
Tgff values below ~27.5 kK; see Vink et al. 2000), we restrict 
our Tgff range to a minimum value of 30 kK. We find that for 
all Fg values M is not a strong function of temperature. In terms 
of the terminal velocity dependence, Fig.[8]shows a rather steep 
dependence on temperature, with Voo dropping by a factor two. 
This is merely a reflection of the escape velocity dropping by a 
similar factor of two over the temperature range under consider- 
ation. 

4.5. Effect of the helium abundance on high Fg models 

To establish the existence of a potential helium dependence on 
M, we computed additional models across the entire Fg region, 
setting the hydrogen abundance to zero and increasing the he- 
lium abundance accordingly. The results are listed in Table |2] 
and shown in Fig. |9] The mass-loss rates are similar to those 
of H-rich models for objects with the same Fg (see Table [TJ. 
This is not too surprising given that the indirect effects of dif- 
ferent continuum energy distributions for H-rich versus H-poor 
are rather subtle (Vink & de Koter 2002). However, when the 



snoo h I I I - 

60 50 40 30 20 

Teff(kK) 

Fig. 8. The predicted terminal velocities versus effective temper- 
atures for several values of Fg, with Fg equal to 0.90 (model 
29; open square), 0.87 (model 28;open triangle), 0.83 (model 
25;open diamond), and 0.58 (model 23; asterisk), respectively. 

He-rich results are plotted in Fig. |9]they lie above the H-rich 
models. For equal luminosity and Fg the masses of the He-rich 
models are lower since Fg is a function of the chemical compo- 
sition through cTg (see Eq. [TJ; cTg is lower for He-rich models, 
therefore the mass must be lowered to keep Fg constant. Similar 
to the H-rich models, there appears to be an upturn in the mass 
loss vs. Fg dependence for models at about Fg ~ 0.7. 

With respect to the terminal velocity and y6-dependence, we 
do not find any significant differences between H-rich and He- 
rich models (see Table [1] versus Table|2). 

5. Spectral morphology: the characteristic He 4686 
Angstrom line 

In the previous section, we provided evidence for a natural tran- 
sition in the mass-loss-Fg exponent, as well as in the velocity 
parameter /? and wind-efficiency from moderate Fg "optically 
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model 
number 


[Mo] 


Wo] 


logL 




[km s"'] 


logM 
[Moyr'] 


/3 


2He 


60 


35.8 


6.0 


0.43 


4552 


-5.36 


1.12 


5He 


40 


23.1 


6.08 


0.80 


3141 


-4.86 


1.42 


lOHe 


85 


49.0 


6.25 


0.55 


4567 


-5.04 


1.21 


14He 


85 


49.0 


6.39 


0.77 


4144 


-4.69 


1.40 


24He 


120 


69.5 


6.50 


0.70 


4800 


-4.70 


1.32 


26He 


180 


104.1 


6.76 


0.85 


4759 


-4.28 


1.53 


29He 


275 


159.5 


6.97 


0.90 


4958 


-3.99 


1.70 


30He 


300 


175.0 


7.03 


0.95 


4934 


-3.88 


1.75 



Table 2. Helium-enriched mass-loss predictions. All parameters that are not listed are the same as in Table[Tl and masses have been 
lowered to keep Tg fixed. 
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log(Gamma) Wavelength (Angstrom) 



Fig. 9. Same as Fig.|5]but for He models from Table |2]added as 
open straight squares. The He models form a straight line above 
the H-rich relationship. 



thin wind" cases to "optically thick wind" cases for objects that 
find themselves above Fe ^ 0.7. We have inspected our models 
and confirmed that for Fe ^ 0.7 the sonic velocity is reached out- 
side the photosphere, whilst the stars form a pseudo-photosphere 
for > 0.7. 

We expect that the occurrence of a pseudo-photosphere has a 
consequence for the spectral morphology of the stars in question. 
We might suspect that the transition Fe = 0.7 is the point where 
the spectral morphology of normal O stars changes from the 
common O and Of-types into a WN-type spectrum. The spec- 
tral sequence involving the Of/WN stars has a long history (e.g. 
Conti 1976, Walbom et al. 1992, de Koter et al. 1997, Crowther 
& Dessart 1998) but it still has to be placed into a theoretical con- 
text. Figure [To] shows a sequence for the predicted He ii 4686 A 
lines for three gradually increasing values of Fg: 0.70 (model 24), 
0.84 (model 15), and 0.93 (model 20), respectively. These mod- 
els have been selected to be objects with a constant luminosity of 
logiL/Lo) - 6.5, and we simply lowered the mass from I2OM0 
to IOOMq, to 9OM0. It is insightful to note that, although the 
first spectrum below the transition Fe already shows some emis- 
sion - characteristic of Of stars - the line-flux is rather modest 
in comparison to what is found for the next two cases with Fg 
values exceeding the critical value of 0.7. These objects show 
very strong and broad He 11 4686A emission lines that are more 
characteristic of Of/WN or "slash" stars, progressing towards the 
Wolf-Rayet stars of the nitrogen sequence (WN). 



Fig. 10. The predicted normalized He n /i4686 flux versus wave- 
length for three values of Fg, with from top to bottom Fg equal to 
0.93 (model 20; 9OM0; EW ~-20A), 0.84 (model 15; IOOMq; 
EW —7 A), 0.70 (model 24, 120Mo; EW ~-2A), respectively. 
Wind clumping has not been taken into account here. 



These models thus indicate that the observed spectral transi- 
tion from Of to WN corresponds to a transition from relatively 
low Fg to high Fg values (and larger /?) for WN stars. This as- 
sertion is based not only on the higher predicted mass-loss rates 
themselves, but also on the finding that the mass-loss behaviour 
(as a function of Fg ) changes at Fg = 0.7. We note that the in- 
creasing Hen 4686A equivalent width (EW) amounts to EW val- 
ues of -2, -7, -20A, respectively (for these unclumped models). 

6. Discussion 

6.1. Comparison with empirical mass-loss rates 

Comparing our new mass-loss predictions against observed 
mass-loss rates is a non-trivial undertaking, as high Fg objects 
are scarce. The largest sample of such potentially high Fg objects 
that involves state-of-the-art modelling analysis is probably that 
of the Arches cluster by Martins et al. (2008). They provided 
stellar and wind properties (accounting for wind clumping) of 
28 of its brightest members from /T-band spectroscopy. Roughly 
half of their sample comprises 04-06 supergiants whilst the 
other half includes H-rich WN7-9 stars. 

It is not possible to quote direct mass-loss predictions, be- 
cause the Martins et al. analysis did not yield object masses. 
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However, on the basis of the high stellar luminosities, with 
logiL/Lo) up to 6.3, these objects were suggested to be consis- 
tent with initial masses of up to ~120Mq. For the 04-6 super- 
giant population, luminosity values are in the range logiL/Lo) - 
5.75-6.05, consistent with initial masses M ^ 55 - 95Mq. For 
this mass and luminosity range, Fg - 0.2 - comfortably within 
our low-Fe regime. Assuming the current mass of these objects 
is about the same as their initial mass, our mass-loss formula 
yields values of log M ^ -6.1, which is in reasonable agreement 
with the lower end of the Martins et al. mass-loss rates for their 
04-06 I objects. 

The second group of Martins et al. objects comprise the 
WN7-9 objects. If we again assume their current masses can di- 
rectly be inferred from the observed luminosities, we find Fg - 

0. 4 and mass-loss rates log M ^ -5.3. Even if the helium abun- 
dances of these objects are increased, these properties would not 
result in a pronounced emission profile of the He ii /14686 line, 

1. e. a profile shape that is typical of late-WN stars. For this to 
happen the mass-loss rates need to be higher by at least a fac- 
tor of a few. This seems to require high F. In the framework 
of our models, this could be achieved by lowering the mass. 
However as the non-electron contribution in our models is not 
self-consistently treated in the high Fg regime, we refrain from 
providing quantitative assessments. 

As the 04-06 supergiants from the Martins et al. (2008) 
analysis have effective temperatures in the range 32-40 kK, we 
expect Ilea to fall in the range 2500-3500 kms"' (see Fig. [8]), 
which reasonably agrees with the upper end of the Arches 04- 
06 supergiant stars. However, for the late WN stars, the termi- 
nal velocities presented by Martins et al. drop to values as low 
as 800- 1 600 km s" ' , which is a factor 2-4 lower than we predict. 
We identify a number of possible reasons for this discrepancy. 
One option could be that the K-band spectral fits of Martins et 
al. (2008) yield terminal velocities that are too low (note that 
no ultra violet P Cygni blue edges are available for these ob- 
scured objects), but a more plausible reason is that as a result of 
our modelling assumptions (no rotation, smooth winds, etc.), we 
overpredict the wind terminal velocity. As the overprediction in 
terminal velocity is a factor 2-4, this would naturally translate in 
a mass-loss underestimate of a factor 3-5. 

Although the effects of rotation on our mass-loss predic- 
tions are beyond the scope of this investigation, rotation may 
become relevant once one wishes to compare the non-rotating 
predictions to observations of objects that might rotate at a rel- 
evant rate. In particular, given that Grafener et al. (2011) dis- 
cuss the possibility that the WNh stars in the Arches cluster may 
evolve close to being chemically homogeneous. To first order, 
one might expect the effect of rotation to lower the effective 
gravity; i.e. one could replace the mass by the effective mass 

M(l -F)(l -Q2), in which case one would anticipate M to scale 
^ , i_ 

as M(Q) oc M(0)(1 - D. ) »' . In a similar vein, one might ex- 
pect the terminal velocity to scale with the escape velocity, i.e. 

Voo Voo(O) V(l -f^^fl (see e.g. Gayley 2000; Puis et al. 2008 
and references therein). In Sect. 14.51 we computed some test 
models in which we lowered the stellar mass, hence the effective 
gravity, so as to keep Fg constant (in order to study a potential 
helium dependence on M). These models indeed showed higher 



^ The issue of rotation in radiation-driven wind is highly complex. 
All existing studies (e.g. Friend & Abbott 1986, Bjorkman & Cassinelli 
1993, Petrenz & Puis 2000, Pelupessy et al. 2000, Cure & Rial 2004, 
Madura et al. 2007) have only treated part of the problem, but not tack- 
led the combined multi-dimensional, high Q (fl S 0.75), and high F 
aspects. 
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mass-loss rates for lower effective gravity, but can we quantify 
the effective-gravity effect? 

If we were to attribute the offset between the observed and 
predicted wind-terminal velocity of a factor ~2 to stellar rota- 
tion, we should have an Q of ~0.85 to "match" theory to obser- 
vations. In that case, the mass-loss rate would be enhanced by a 
factor of 2.5 - 4 for a' values in the range 0.5-0.6. Nevertheless, 
this could significantly change for lower temperatures. These 
high Q values of the order of 0.85 seem to be rather high when 
we consider recent stellar evolution models (Brott et al. 2011; 
Friedrich et al. in prep.), but they are not out of line with respect 
to spectral observations of LB Vs (Groh et al. 2006), objects that 
are presumably in close proximity to the Eddington limit. 

In all, we interpret our high wind terminal velocities as a sign 
that some physics is missing in our high F models. Therefore, we 
refrain from using our dynamically determined mass-loss rates 
in a quantitative way at the optically-thick high F end. For nor- 
mal O-type stars, for low and moderate Fg values, we achieve 
much better agreement between observed and predicted wind 
terminal velocities (Muijres et al. 2011b), and for this regime 
we have a much higher confidence in the accuracy of the abso- 
lute mass-loss rates, as long as O-type winds are not extremely 
porous, in which case mass-loss rates could drop significantly 
(Muijres et al. 2011a). 

6.2. Comparison to other models 

6.2.1 . Comparison to CAK and other O-type star mass-loss 
models 

We now wish to compare our results with previous model predic- 
tions. In this paper, we have investigated the mass-loss behaviour 
at high Fg for an extensive grid of models, and we revealed the 
existence of two mass-loss regimes. Moreover, we have found 
that mass-loss rates are dependent on both Fg and stellar mass 
(or stellar luminosity) and that the shape of these dependencies 
is well-described by a power law. 

We remind the reader of classical CAK theory, where the 
mass-loss rate is found to be proportional to 

M oc L (4) 

where a is a force multiplier parameter expressing the impor- 
tance of optically thin lines to the total ensemble of lines. It is 
generally found that a is ~2/3 for galactic O-type stars (Puis et 
al. 2008) and assumed to be constant throughout the atmosphere. 
In reality, however, a is depth-dependent (Vink 2000, Kudritzki 
2002, Grafener & Hamann 2005, Muijres et al. 201 lb), which is 
captured better by an alternative representation of the line accel- 
eration (Miiller & Vink 2008). Nevertheless, the classical CAK 
formalism - as described by Eq. |4] - already shows a depen- 
dence on both L and Fg, and one could rewrite this mass-loss 
dependence as a function of M and Fg, using a mass-luminosity 
relation. 

In the standard Vink et al. (2000) mass-loss parametrization 
M oc l2.2 ^-1.3 (y^/u^^j-i-z^ which can be reorganized to M oc 
^1.2 rpj^jg jypg rnass-loss parametrization that is 

currently employed in modem evolutionary computations (see 
e.g Meynet & Maeder 2003, Palacios et al 2005, Limongi & 
Chieffi 2006, Eldridge & Vink 2006, Vink et al. 2010, and Brott 
etal. 2011). 

In this paper, we obtain much steeper M vs. F dependencies, 
in agreement with our previous models for constant-luminosity 
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LBVs (Vink & de Koter 2002, Smith et al. 2004). For the "low" 
Te range considered here we obtained "modest" dependencies 
of M oc l"-^ Fe''^ , but for the "high" Fg regime, we found M 
oc L"-^^ Fe"^^^ which involves a much steeper dependence on Fe 
than any CAK-type mass-loss relationship provides. We note 
that such a steep dependence on the Eddington limit agrees with 
radiation-driven wind models of Vink (2006) and Grafener & 
Hamann (2008), whilst Grafener et al. (in prep.) also provide 
empirical evidence for such a strong mass-loss dependence on 
the Eddington parameter. 

We emphasize that, through the use of the Vink et al. (2000) 
theoretical mass-loss recipe, most current stellar models already 
include the effect of positive mass-loss feedback (contrary to 
recent claims by Smith (& Conti 2008). This effect describes 
how the mass-loss rate increases with the Eddington parameter. 
However, as we here obtain much steeper M vs. F dependencies, 
it is likely that the mass-loss feedback effect that is currently 
employed in the stellar evolution models may not be sufficient 
for certain areas of the Hertzsprung-Russell diagram. We thus 
concur with the notion of Smith & Conti (2008) that new stellar 
evolution computations that take this effect into account properly 
are desirable 

6.2.2. Comparison to alternative Wolf-Rayet mass-loss 
models 

We also compare our models to the optically thick wind models 
for Wolf-Rayet stars, such as the ciitical-point analysis by Nugis 
& Lamers (2002) and the hydrodynamical model atmosphere 
analysis of Grafener & Hamann (2008). As there is a signifi- 
cant qualitative difference between our Monte Carlo approach 
and the optically thick wind approaches, a meaningful quantita- 
tive comparison is a non-trivial undertaking (see section ISTTT i. 

First, we quantitatively compare our M versus Fg depen- 
dence to the WNL star mass-loss dependence suggested by 
Grafener & Hamann (2008). For the models in our grid at Teff 
= 50 kK, we find very good agreement with the Grafener & 
Hamann (2008) mass-loss rates and also find that the power-law 
slope of our dependence is very similar. However, the onset of 
WR-type behaviour occurs earlier, i.e. for lower Fg, in the mod- 
els by Grafener & Hamann. 

In Sect. [3] we discussed the possibility of such a shift in Fg, 
because the actual Eddington parameter F is expected to be af- 
fected by free-free and bound-free contributions and peaks in the 
iron opacity. By comparison with OPAL opacity tables (Iglesias 
& Rogers 1996), we estimate an increase of F by ~ 20% in the 
region of the sonic point, assuming the location of the sonic point 
remains unaffected. This value corresponds roughly to a ~ 25% 
shift in Fg between our relation and that by Grafener & Hamann 
for typical parameters of Galactic WNL stars {T^s - 45 kK, 
log(L/L0) = 6.3). This could be considered a maximum shift 
as we may slightly overestimate the line force near the sonic 
point by applying the Sobolev approximation (Pauldrach et al. 
1986). However, a change in F affects the atmospheric structure 
and therefore the location of the sonic point, consequently the 
effect on M is hard to establish. 

The Grafener & Hamann mass-loss rates also display a 
strong temperature dependence, with M oc T^^-^ . The actual size 
of the shift in Fg is thus strongly dependent on the specific stellar 
parameters. Our Monte Carlo models suggest a much smoother 
dependence on Tgff (see Fig. |7]i as long as we stay above the 
location of the predicted bi-stability jump, where the mass-loss 
properties jump drastically (Vink et al. 1999, Pauldrach & Puis 
1990). 



What we wish to emphasize is that both modelling ap- 
proaches show an M versus Fg dependence that is much stronger 
than any additional mass or luminosity dependence. Where the 
two distinct mass-loss prescriptions differ is in the treatment of 
wind clumping, as well as the value of F for the onset of WR- 
type mass loss behaviour. The exact location of this transition 
is of paramount importance for the evolution of the most mas- 
sive stars. Ultimately, this should be testable with comparisons 
to observational data when sufficient objects are available in the 
appropriate Fg range. This will be a crucial aim of the VLT 
Flames Tarantula Survey ( Evans et al. 2011, Bestenlehner et 
al. in prep.). 

7. Summary 

We presented mass-loss predictions from Monte Carlo radiative 
transfer models for very massive stars in the mass range 40- 
300Mo and with Eddington factors Fg in the range 0.4-1.0. An 
important outcome is that when winds become optically thick 
the spectral and mass-loss properties change. This transitional 
behaviour can be summarized as follows 

- We find a transition from common O-type stars to more 
extreme Wolf-Rayet behaviour when Fg exceeds a critical 
value. 

- The way in which the mass-loss rate depends on Fg in the 
range 0.4 :S Fg ^ 0.7 is M oc M*° ''^Fg-2, where rates are found 
to be consistent with the standard Vink et al. (2000) mass- 
loss rates. 

- At Fg ^ 0.7 the M dependence shows a "kink"; i.e. the slope 
is steeper for objects closer to the Eddington limit. Here 
the slope becomes M oc M*"'^^F^'^^. This slope agrees with 
WNL models by Grafener & Hamann (2008). 

- When Fg approaches unity, the wind efficiency number rj 
rises in a curved manner to values as high as rj ^2.5. 
Such high T] values are more commensurate with Wolf-Rayet 
winds than with common O stars winds, and these results 
thus confirm a natural extension from common O-type mass 
loss to more extreme WR behaviour. 

- This transitional behaviour is also found in terms of the wind 
acceleration parameter y6, which naturally reaches values as 
high as 1.5 

- The spectral morphology of the He ii line at 4686A changes 
gradually as a function of Fg. This links the spectral sequence 
O-Of-Of/WN-WN to a transition of optically thin to opti- 
cally thick winds. 

- The mass-loss rate is found to be only modestly dependent 
on the eff'ective temperature for the range of 30 to 50 kK. 

- Last but not least, we highlight the fact that for fixed Fg the 
He abundance only has a minor effect on the predicted rate of 
mass loss (cf. Vink & de Koter 2002 for LBVs). This contra- 
dicts how O-type, to LBV-type, to WR-type mass-loss tran- 
sitions are employed in massive star evolution models. We 
thus call for fundamental changes in the way mass loss is in- 
cluded in stellar evolution models for objects in close prox- 
imity to the Eddington limit. 

Acknowledgements. We would like to thank the anonymous referee for provid- 
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